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ABSTRACT 

Femtolensing is a gravitational lensing effect in which the magnification is a 
function not only of the positions and sizes of the source and lens, but also of the 
wavelength of light. Femtolensing is the only known effect of (10"^^ — 10~^^ Mq) 
dark-matter objects and may possibly be detectable in cosmological gamma-ray 
burst spectra. We present a new and efficient algorithm for femtolensing 
calculations in general potentials. The physical-optics results presented here 
differ at low frequencies from the semi-classical approximation, in which the 
flux is attributed to a flnite number of mutually coherent images. At higher 
frequencies, our results agree well with the semi-classical predictions. Applying 
our method to a point-mass lens with external shear, we find complex events 
that have structure at both large and small spectral resolution. In this way, we 
show that femtolensing may be observable for lenses up to 10~^^ solar masses, 
much larger than previously believed. Additionally, we discuss the possibility of 
a search for femtolensing of white dwarfs in the LMC at optical wavelengths. 

Subject headings: dark matter — gamma rays: bursts — gravitational lensing 
— methods: numerical 



- 3 - 



1. Introduction 



The possibility of interference effects in gravitational lensing has been considered by 
several authors ( [Mandzhos 1981| ; [Schneider fc Schmidt-Burgk 1985| ; peguchi fc Watson] 



[I986t [Peterson fc 1^'alk 1991| ; [Gould 1992| ; [Stanek, Paczyhski, fc Goodman 19931 ). Of 



particular interest are diffractive variations in flux with frequency v when the source, lens, 
and observer occupy fixed positions. If the difference in time delay between a pair of 
images satisfies vtd 1, the fringe spacing is Az/ ^ At"^ ^ v. We may call z/At ^ 1 
the "semi-classical" regime, because diffractive phenomena are produced by mutually 
coherent images whose positions, magnifications, and time delays can be determined using 
geometric optics. On the other hand, if z/At <^ 1, regions of the lens plane other than the 
geometric-optics images contribute importantly to the fiux. The semi-classical approach 
then breaks down, and one must use the methods of physical optics. 

If observed, interference effects would reveal dark-matter objects in a mass range to 
which few other tests are sensitive. The characteristic time delay produced by a lens of 
mass M is 

At(M) = 2GM/c3 = i?sch/c, (1) 

where -Rsch is the Schwarzschild radius. Hence the condition z/At ~ 1 is equivalent to 
A ~ -Rsch- Since the Schwarzschild radius of the Sun is ~ 3km, broad-band fringes 
(Az/ ~ v) require decidedly sub-stellar but nevertheless macroscopic lensing objects. In 
particular, Gould (1992) has shown that lens masses M ~ 10^^^ - lO^^^M© ~ 10^^ - lO^^g 
could produce observable fringes in gamma-ray burst spectra at energies E ~ IMeV 
(A ~ lO^^'^cm). Because the angular separation of the images produced by such a lens is 
~ 10^^^ arc sec, Gould has coined the name "femtolensing" for this phenomenon. On the 
other hand, the probability that a randomly placed and cosmologically distant point-like 
source should be lensed is ~ f^iens, where the latter is the mean mass density in lensing 



-4- 



objects expressed as a fraction of the critical density SHq/SttG ( [Press fc Gunn 1973|) . 
Given that ~ 10'^ gamma-ray bursts have been detected to date (CGRO Science Report 157 
1994), a single well-established case of femtolensing would indicate that objects ~ IO^^^Mq 
contribute significantly to the mass density of the universe. 



Even if copious lenses exist in the appropriate mass range, visible fringes can be seen 
only if the (incoherent) source is smaller than the Fresnel length (AD)^/^, where D is 
the distance. For source redshifts of order unity, this translates to -Rsource < lO^^Acmcm. 
An additional constraint requires that in order to have significant magnification, the 
source size must be smaller than the Einstein ring radius. For cosmological distances, 

i^source < 5 X lOHl / {10'^^ Mq) Cm. 

The rapid time variability of gamma-ray bursts (GRB) of 0.2 msec (IBhat et al. 1992| ) 



as well as the cosmological distances of ~ 0.5 Gpc for bright BATSE bursts found from 
Log A^-Log P studies ( [Fenimore et al. 1993| ) suggest, however, that GRB are sufficiently 



compact. There has been some confusion as to whether the appropriate linear source size 
should be taken from 7cAt or 7ct where 7 is the bulk Lorentz factor, c is the speed of 
light. At is the smallest time variation detected, and t is the total event duration. For 
many proposed cosmological scenarios, the appropriate measure is 7cAt, because the 
last-scattering-surface remains at approximately the same radius even though the relativistic 
ejecta may reach quite large distances in the course of the burst. For 7's of 100-300 (e.g. 
b'enimore, Epstein, fc Ho 1991] ), the (linear) size of a GRB last-scattering-surface is of order 



5 X 10® cm. Other models predict emission from patches on a relativistically expanding 
shell which can become extremely large (7ct). In the latter case, it would be difficult to 
observe femtolensing because the interference patterns would differ from patch to patch. 
An observation of femtolensing could distinguish between the two scenarios. 

Stanek, Paczynhski, & Goodman (1993, henceforth SPG) have discussed the possibility 
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that line features in burst spectra may have been produced by femtolensing. Such hues 
have been seen or inferred in GINGA data ( [Murakami et al. 1988| , [Fenimore et al. 1988|) and 



KONUS data ([Mazets et aL 19^ ) and have been attributed to cyclotron absorption. Like 



cyclotron lines, interference fringes would be evenly spaced in photon energy. No convincing 
evidence for lines has yet been seen in the largest homogeneous data set available, the 
BATSE experiment on the Compton Gamma-Ray Observatory, although its capability of 
line detection is lower ( [I'eegarden et. al. 1993| ). 

The femtolensing calculations cited above have considered only the simplest possible 
case, which is an isolated point-mass lens. The simplicity and symmetry of such a lens 
allow the physical optics problem to be solved in terms of known functions (DegucEFS 



Watson 1986| ). In the present paper, we present an efficient physical-optics method for 



computing frequency-space fringes produced by general lenses. We assume that the lensing 
mass distribution is confined to a layer thin compared to the observer-source distance 
(single-screen approximation). We also neglect time dependence of the lensing geometry, 
which is permissible if the time-delay difference between any pair of images changes by less 
than z/~^ during an observation. The latter assumption is probably justified for femtolensing 
of gamma-ray bursts ( [Gould 1992| ). 



For definiteness, and because it is the simplest lens not yet treated in physical optics, 
we apply our methods to a single point mass with external shear. The computational 
approach taken, however, would apply equally well to an arbitrary surface density of lensing 
mass. Computational savings are achieved mainly by taking advantage of the achromaticity 
of gravitational lensing: that is, the time delays and excess optical path lengths are 
independent of frequency. 

The plan of our paper is as follows. The physical-optics problem is posed in §11. We 
show how the scalar diffraction amplitude at the observer can be determined as a function 
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of frequency by first calculating its Fourier transform, which is a function of time delay. An 
efficient numerical procedure for finding the latter function is developed in terms of contour 
integrals on the lens plane. The role of the geometric-optics images and the correspondence 
with the semi-classical approximation is explained. In §111 we describe certain details of our 
numerical implementation of the method and show that in the case of an isolated-point-mass 
lens, our results are consistent with those already obtained by [Deguchi fc Watson 198(j| and 
SPG for the isolated-point-mass lens. In §IV we present results for the more complex case of 
a point mass with external shear, which can produce up to four images. Finally, §V briefiy 
summarizes the main points of our approach and prospects for observing femtolensing. 



2. Formalism 



The quantity of interest for femtolensing is the observed magnification as a function 
of frequency of the lensed source relative to the intensity in absence of a lens. In scalar 
diffraction theory, for the case of a thin screen, this function can be written as complex 
square of the amphtude 



oo roo 



(2) 



oo J — oo 



where u is the photon frequency, is a normalization to unit flux in absence of a lens 
that varies slowly with u, x and y are coordinates in the lens plane, and r is the time delay 
function (e.g., [Blandford fc Narayan 198^) : 



r(f) 



1 + Zi 



dn 



J dx'dy'T,(f') In |r — 'P' 



(3) 



_2doidis 

In the latter formula, zi is the redshift of the lens; dos, doi, and dig are angular diameter 
distances; = {xs,ys) is the point where a direct line from observer to source would meet 
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the lens plane in the absence of the lens; and S(r^ is the mass per unit area in the lens 
plane. 

Although our methods are general, we illustrate them by application to a point mass 
with sub-critical shear. If the observer and lens are on the z-axis and the source is slightly 
off axis, the time delay function for this case is, in normalized units (see Appendix A), 



T{x,y, iJ,,(j),9) = (1 + + (1 - - 26'[xcos(0) +?/sin(0)] - ln(a;^ + y'^), (4) 



where /x is the shear, is the angle between shear direction and that of the displacement of 
the source from the axis, and 6 is the angle between the source and z-axis as measured by 
the observer in units of the angle subtended by the Einstein ring. 

Direct calculation of \I^(ci;) is a three dimensional problem in x,y, and u. We 
demonstrate here a method to reduce the problem to two dimensions with use of Fourier 
transforms and contour integration. Dividing Eq. |^ by C^^ and taking the Fourier transform, 
we define 



which can be viewed as a virtual pulse shape in time. After substitution from Equation 0, 
this equation reduces to 



Therefore, the contribution to ^{t) comes from curves of constant time delay; an example 



curves defined by T{x,y) = t and T{x,y) = t + dt, and the distance between them is 




(5) 




(6) 



of such contours is shown in Figure |I[ 



Eq. ^ can be evaluated as a contour integral. Since \E'(t)dt is the area between the 



dt/||Vr 




(7) 
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where ds is the arc length along the curve defined below. 



There will in general be more than one of these contours, though it is straightforward 
to find them all as is described in the next section. The sum of the integrals (0) over all 
contours {7^} can be expressed as 



^(t) = ^ I / du\ 



if the arc parameter u is defined through the differential equations 

dx dr ^ dy dr 
du dy du dx'' 



because 



ds 
du 



^(1)^(1)-"-"- 



(9) 



(10) 



Notice that equations (P) are nonsingular near critical points if r(r) is smooth. 



We find '^{t) is a smooth function except where t is equal to a critical value of the time 
delay function so that ||Vr|| = (an example of \l/(t) is shown in Figure It is at these 
critical valued time delays that images form. These singularities are calculated in Appendix 
B. For minima or maxima, there is a discontinuity in '^{t) at the critical time, tj, of 



'd^Td^T 


( VI 


dx"^ dy"^ 


\dxdy J 



lim{^ {ti + e) - ^(tj - e) = ±2n 



In the case of a saddle point, \l/(t) diverges logarithmically near tj, so that 



'11^ 



*(t) 



-2 In It - t, 



d^T d' 



+ non-singular part 



(12) 



\ dxdy J dx"^ dy"^ 

Using a procedure described in the next section, the singularities can be dealt with and the 
inverse Fourier transform can be calculated, yielding the amplitude 



^'(u;) = Cuj exp(?cjr)^(t). 



(13) 
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3. Numerical Implementation 



Except in the simplest circumstances, tlie metliod of §2 must be applied numerically. 
For a point mass with shear, our algorithm computes the femtolensed energy spectra in 
ten to twenty seconds on a Sparc 10. Computer code which implements the algorithm 
is available upon request. For future applications, we note that the algorithm grows 
quadratically with the number of lensing masses. The computation may be considered in 
three parts: location of the images, line integration on constant time delay contours, and 
calculation of the inverse Fourier transform of \l'(t) to yield \l'(a;). 

In the case of a point mass lens with shear, the location of the images can be found 
analytically, so there is essentially no computational burden. However, one could envision 
more complex applications such as lensing in a dense field with many lenses and images 
where it would be necessary to utilize some image finding scheme. An efficient image finding 
algorithm has already been developed for microlensing calculations ( |Witt 19^31) . Note that 
the femtolensing algorithm increases approximately linearly with the number of images, so 
these complicated scenarios may be tractable. We find that the line integration around 
contours defined in Eqs. p^-jTOl is the most computationally demanding part of the program. 



The method by which the contours are found warrants a brief discussion. An example 
of these contours are shown, for example, in Figure |I[ The algorithm starts with a contour 
very close to the minimum time delay critical point. From there, it moves to larger time 
delay contours by taking small steps perpendicular to the contour. This method works as 
long as the number or topology of the contours does not change, but this happens only at 
the images. In other words, the contours are altered only when moving from below the time 
delay of an image to above it. When the image is a minimum in the time-delay surface, a 
new contour is formed around the minimum. When the contour is a saddle point, a contour 
either splits into two, or two contours merge. One can distinguish between these cases by 
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monitoring the line integral at contours near each time delay of an image. Maxima are not 
encountered in the case of a point mass with shear. 

Using the line integrals we determine \E'(t) for a range of values of r, and the 
computations for one geometry are shown in Figure |^. We find generally that one needs 
approximately 1000 points taken uniformly between the minimum time delay and the 
point at which reaches about one percent of its asymptotic value. It is impossible 
to reconstruct the entire detail of the ^{t) curve with such measurements, because for 
instance, near the time delay of a saddle point image there is a logarithmic divergence as 
discussed in the previous section and Appendix B. We account for this loss of detail near 
the critical time delays by subtracting the analytic forms of the singularities. (Eqs. |TT], |T^) 
and dealing with them separately: 

^'{t) = X! I r '^^1 ~ singular part (14) 

k •'"Ik 

We then perform a discrete Fourier transform (DFT) of while paying close attention 

to edge effects. Finally, we calculate the exact Fourier transforms of the divergences and 
discontinuities and add them to the DFT to obtain '^{u). 

As a test of this procedure, the calculation is carried out in the simplified case of a 
point mass without shear and is shown in Figure ^ (top panel). This calculation agrees with 
an earlier calculation of the same result which utilized independent methods (SPG) to the 
resolution of the figures presented in this paper. 

4. Results 

Using the new algorithm, we investigated the interference produced by a variety of 
lensing geometries for a point source and lens with shear (extended sources could be 
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calculated by repeatedly applying the algorithm for source points chosen, for instance, by 
Gaussian quadrature weighings). Figures ^ and ^ show the results of calculations for various 
source/lens configurations depicted in Figure ^. The frequency, tu, is given in dimensionless 
units: 



where Rg is the Schwarzschild radius and A is the wavelength. For example, unity 
corresponds to photon energy of 1 keV and a mass of ~ 5 x 10~^^Mq at zl = 0.5. A 
sinusoidal pattern is produced which has a frequency inversely proportional to the difference 
in time delay for the two images. The top panel of Figure |^ shows the interference pattern 
for a point mass with no shear for which the semi-classical approximation is an extremely 
good approximation. 

Using our algorithm, external shear can be introduced into the calculation. We find 
that when the source approaches the caustic, the physical optics calculation differs from 
the semi-classical approximation. For example, the bottom panel of Figure ^ shows a two 
image interference pattern for which the first few interference nodes have largely different 
amplitudes. At high frequencies, the semi-classical approximation approaches the physical 
optics solution. 

The interference patterns become much more complex when the source lies inside 
the center caustic region. There, four images form, and the interference patterns become 
chaotic. Figure |^ shows the magnification as a function of frequency for two such cases. 
The main characteristics to note are that the interference patterns can become quite 
complex with oscillations with both long and short periods (the bottom panel shows both 
the physical optics magnification and a smoothed magnification). In short, this means 
that femtolensing can be observed over a wider range of frequency space, and therefore, a 
wider range of masses than originally believed. Figure |^ shows that the typical interference 




(15) 
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patterns, corresponding to inverse time delays, can be hundreds of times larger than the 
characteristic time delay (see Eq. |l]), so the lens masses can be hundreds of times larger than 
those discussed by Gould and SPG. This short period interference discussed in previous 
papers can be seen in the top panel of ^ The long period interference occurs when two of 
the four images have very close time delays, which is a common result with four images. 
Appendix C shows that in these circumstances the separation of the long-period fringes 
increases as the cube of the total magnification, and the source-size limit is somewhat more 
severe than the Fresnel length. 

The complex nature of the magnification in Fig. 5 can be understood to some 
extent in the semi-classical regime, because there one finds that the complex amplitude of 
magnification, "^{w), is composed of a set of sinusoidal interference fringes from each pair of 
images. For four images, there are six pairs which each produce sinusoidal interference. The 
resulting interference pattern which is the squared sum of these sine waves will in general 
be very complex. The physical optics calculation is qualitatively similar, but yields different 
magnifications and different specific structures at low frequencies. According to the result 
(^) in Appendix C, the first few long-period fringes are least sensitive to the angular size 
of the source. 

5. Discussion 

We find that even in the simplified case of a point mass with external shear, very 
complex interference patterns can be formed. In reality, however, the patterns should be 
even more complex. In addition to external shear, one should account for local shear from 
neighbors for femtolensing matter in galactic halos. Complex, many image geometries 
will result, and as found in the case of a point mass with external shear, images with 
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smaller-than-characteristic difference in time delay (see Eq. |l|) will cause surprisingly long 
period interference patterns (e.g. Figure ^ which allow for detection of larger masses. 
For an external shear of 0.1, we expect these complex interference patterns to occur 3% 
of the time that a source is found inside the Einstein radius. In general, for a halo mass 
distribution similar to a singular isothermal sphere, the shear would have a scale length 
of 5 kpc, so that between 10 and 50 kpc, the shear would be 0.5-0.1, so these complex 
patterns would occur much more often. Furthermore, there would be a magnification bias 
(see Appendix C) towards these complex events, which could likely increase the observed 
fraction by an order of magnitude. 

Cosmological gamma-ray bursts are considered the best candidates to show femto- 
lensing because of their extremely small angular size. If the line of sight to a gamma-ray 
burst passed through a galaxy, and if the dark halo mass are composed of ~ 10^^^ — IO^^^Mq 
(we increase the range of masses by a factor of 100 as a result of the long period events 
discussed above) then one could expect to see, in addition to macro or microlensing, 
femtolensing effects as well. The femtolensing could easily mimic other emission or 
absorption line processes. We suggest, then, that if a spectral absorption or emission feature 
is definitively observed in a gamma-ray burst due to femtolensing, that it is likely that the 
burst would be macrolensed as well. Additionally, macrolensing would produce multiple 
images (bursts) which would likely be femtolensed and therefore increase the number of 
expected femtolensed events. 

Although the source size requirements exclude the possibility of ever detecting 
femtolensing in main sequence stars in nearby galaxies, in the future, it may be possible 
to detect femtolensing of white dwarfs. (We are grateful to B. Paczyhski for pointing this 
out.) In order for the fringes to have separation of 0.2 eV (a typical optical bandwidth), the 
lensing masses would be in the range M ~ 10^^ — IO^^^Mq for short to characteristic time 
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delays, respectively. In the LMC, the maximum magnification possible for such events as 
determined by the ratio of the Einstein ring to source size would be ~ 2 — 20 where the larger 
masses cause the higher magnification. For a white-dwarf source in the LMC observed at 1 
eV, the requirement that the source be smaller than the Fresnel length, (A-D)^/^ = 3 x 10^ 
is met {D = Dlmc/'^ = 25kpc). These events would be quite short (< 1 minute), due 
to relative velocities and a small Einstein ring, relative to current microlensing events 
(e.g. [Paczyhski 19861). However, it is not yet possible to monitor white dwarves in the 



LMC, so one would have to adopt an observing strategy for transients — possibly in the 
ultraviolet — similar to that of supernovae searches. In contrast to current microlensing 
studies, detection techniques would be based on anti-correlation of the flux variation in 
different wavelength bands. 



This work was supported in part by the David & Lucille Packard Foundation and 
NASA grant NAG5-1901. It is a pleasure to acknowledge B. Paczyhski for suggesting 
the white dwarf lensing scenario as well as E.E. Fenimore, K.Z. Stanek, and H. Witt for 
discussions. 
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6. Appendix A: Normalizations 



In this appendix, we describe tlie normalized units used in the formulae above. We 
begin with the time-delay formula for a point mass with shear (e.g. piandford fc Narayan 
19921) : 



~, X l+zi 



(16) 

where zi is the lens redshift, x,y are points at which the ray intersects the lens plane, Xs-,ys 
give the location of the source in the source plane, -Rsch is the Schwarzschild radius, and 
D = didis/dg where d\,ds,d\s are angular diameter distances. 

Throughout the derivation, additive constants will be ignored. The form for \l/(co') is 

_ poo peso 

■^{u)=Cco / dxdyexp{iCuf). (17) 



— oo J — oo 



Then, we define the following dimensionless variables: 

fc 



time delay: r 



(1 + zi)Rsch 



u;{l + ZL)Rsch 2tt{1 + ZL)Rsch 

frequency: oo = = (19) 

c A 

(x, y) 

position (image plane): {x,y) = . ' = (20) 

V 2-Ditsch 

shear: /i = 2DRschfi] < fi < 1. (21) 

With these definitions Eq. ^ becomes Eq. ^ and C,^ which normalizes the magnification 
to a unit flux, can then be shown to be, 

a = ^ (22) 

TT 



Equation IS can be rewritten after substituting the dimensionless variables as 



t{x, y, fi, Xs, ys) = (1 + f^)x^ + (1 - - ^ '^DjxsX + y^ _ ^^^^^ ^ ^^^^ ^^g) 
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Finally, we choose an impact parameter 



e 




(24) 



where i?E is the Einstein radius: 



-Re = 




(25) 



We also pick cf) so that we recover Eq. ^ 



T{x,y,^,4>,9) = (1 + /i)x^ + (1 - fi)y'^ - 26'[xcos(0) +?/sin(0)] - ln(x^ + ?/^). (26) 



Herein, we discuss the behavior of \l'(t) near critical points (images), so that the 
divergent portions of the function can be subtracted to yield a smooth function that can be 
computationally sampled and transformed. When the frequency is larger than the minimum 
separation of image time delays, ie. uj ^ ki ~ '^jl ^lU i, j, it is in the semi-classical region 
and these singularities dominate "^{u). In any regime, the singular behavior dominates ^(t) 
near the time delays of the images. 

We consider the two cases of a minimum/maximum and a saddle point: 



7. Appendix B: Semi-Classical Approximation 



7.1. Minimum/Maximum 



Near a critical point at the origin, a contour is given by 




(27) 
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where tcrit is the time delay at the critical point x = y = 0. Without loss of generality we 
consider the case of a minimum only. The contribution of the singularity to "^{t) is found 
by the integral (see eqs. 



We write: 



dv? 



dx dr 
du\ where — = —7—, 
ou oy 



d dr d V 

X = A cos(-^ + 0) 
ab 



dy dr 
du dx 



X 



so that 



and 



Further, from eq. 27 



du = 2'Kah. 



a^6^ dx^ dy'^ 
The singularity is such therefore, that 



^dxdy 



? = detH. 



if t < tmir 

271 /V detH if t > 



(28) 

(29) 
(30) 
(31) 

(32) 



(33) 



This term can be removed from the slowly varying part of \E'(t) and Fourier transformed to 
yield the semi-classical contribution as 



A-^{uj) = / dtexp{iujt)A'^ 



VdetH 



(34) 



7.2. Saddle Point 



In a similar manner as above, the singularity of a saddle point may be calculated. 
Without loss of generality we consider this orientation of saddle point: 



T 



^2^~ 2^' 



(35) 
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By analogy with the last derivation, we have 

X = Acosh.{u/ab) (36) 

In this case the contours are hyperbolic rather than elliptical and do not close locally. To 
compute the contribution from the neighborhood of the critical point, we consider that part 
of each contour for which |x| < Xo, where Xo is fixed as r — tcrit- Finding Uo such that 
x{Uo) = Xo, and taking \uo/ab\ ^ 1, we have 

Uo = a6cosh^"'^(— ^) -ablog{-^) = — afolog \t — tcriti + constant, (37) 

where the constant depends on the choice of Xo- The contour extends from —Uo to Uo on 
each of 2 branches of the hyperbola, so 

A^{t) = -2abhg \t - t„it| = (38) 
^' ' V- det H ^ ' 

When Fourier transformed, \I^(t) yields the semi-classical contribution, 

Avl/(^) = cj dtexp(z^t)A^ = 2^^^P(^^^cnt) ^3g^ 
J V — det H 



8. Appendix C: Physical Optics of Merging Images 

When the flux is dominated by two bright images near a critical line, the time-delay 
surface can be approximated locally by 

x^ 

r{x,y]Xs,ys) = ^ + - XsX - y^y. (40) 

Here {xs,ys) is the undeflected source position in the lens plane. We use the dimensionless 
units of Appendix A, but the local coordinates in equation (^Dj) are centered on the critical 
line {x = 0) rather than the deflecting mass. The parameters a and b are assumed to be 
~ 1, but their precise values depend on the global geometry of the time-delay surface. 
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For Xs > 0, geometric optics yields two images at x = ±a"^/^xy^ with net flux 

Fgeo(x,)=2a3/262a;;i/2. (41) 

On the other hand, the amplitude can be evaluated directly: 

F{x„u) = \^{u-x,,y,)\'=A7ruj'W\Ai{-aXsU;'/')\'. (42) 
The Airy function quickly approaches the asymptotic approximation [e.g., [Abramowitz fc 



Stegun 1972 



Ai(-0 ^ vr-V^^-^/Sin + ^) ife>l, (43) 

so that the flux (^) displays fringes with spacing 

Au;f=^{ax,r'/'; (44) 

in fact, even the position of first zero in the flux is accurately predicted by the asymptotic 
form (|3), which is the semiclassical prediction. Both geometric optics and the semiclassical 



theory predict a divergent flux as 0. But since Ai(0) = 3-2/3/r(2/3), the flux at the 

caustic is finite and scales as cu^^^. 

Comparing (^) and (^^, we see that the fringe spacing scales as the cube of the total 
magnification. This accounts for the long-period fringes that we often see in our numerical 
spectra. Because of magnification bias, observed fringes should typically be much more 
widely spaced than the naive scaling Auj ~ c/Rsch assumed in earlier work. 

A source of finite angular size can be regarded as an ensemble of incoherent points. If 
fringes are to be visible, the phase of the Airy function in ( ^2]) must vary by less than the 



fringe spacing (^^ across the source. Therefore the angular size must satisfy 

«.«9.^f^)^-n (45) 



where 6e = y'2DRsch/di is the angular radius of the Einstein ring. For images close to a 
critical line, this limit supersedes the angular-size limit associated with the Fresnel length, 
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which scales as u"^!'^. To achieve fringes of a given physical frequency Co and spacing Ao; 
with lenses whose Schwarzschild radius is much larger than the wavelength, the limit on the 
source size scales as -Rsch*^- 
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Fig. 1. — This figure shows the contours of equal time delay in the image plane that are 
used to calculate ^'(t) (cf. fig. The empty circle in the middle contains the lens and 
therefore a logarithmic spike in time delay. Images form at critical points on this surface. 
Two minima form contours at low time delay on the top and bottom right of the depicted 
region. The contours merge at a saddle point on the right. This contour later splits to form 
inner and outer circular contours. 
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Fig. 2. — Figure § shows ^'(t), which yields the magnification amphtude, "^{u) when Fourier 
transformed. In this four image geometry ^{t) has four singularities (though the first occurs 
at minimum time delay). Two minima at low time delay cause discontinuities. Two saddle 
points at higher time delay form logarithmic spikes. The inset shows \E'(t) in a case with only 
two images as the source approaches a caustic. The function resembles that of the four image 
case. This transition region is unique to the wave nature of the physical optics calculation 
and contrasts with the semi-classical approximation which has abrupt qualitative changes 
between two and four images. 
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Fig. 3. — Shown in the image plane are the image positions (crosses), critical curves 
(circle and ovoids), projected source positions (large dots), projected caustics (small dot 
and diamonds) in dimensionless units for a variety of geometries. Fig. ^ is a point mass 
without shear, so the caustic is the Einstein ring. Fig. P3,c,d have shears of 0.1, 0.3, and 
0.1. 
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Fig. 4. — Figure § shows the magnification as a function of dimensionless frequency 

for the geometries in Fig. 0a,b (top and bottom respectively). The sohd hne is the physical 
optics calculation calculated with our algorithm. The dashed line is the semi-classical, or 
two-ray, approximation discussed in appendix B. When the source approaches a caustic, the 
energy spectrum cannot be accurately calculated with the semi-classical approximation at 
lower frequencies. 
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Fig. 5. — This figure shows magnification as a function of dimensionless frequency for the 
geometries in Fig. |^,d (top and bottom, respectively). Complex beat frequencies are created 
between the time delays of the four images so that the semi-classical approximation (dashed 
line in top panel) is inaccurate. The beat frequencies create features on a much broader 
energy range as exemplified by the bottom panel in which a smoothed magnification curve 
is superposed on the physical optics calculation. 



